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1. Introduction 

1.1. Antecedents. In virtually all known applications of the Ricci flow, it is valu- 
able to have a good understanding of singularity formation. Heuristically, there 
are at least three reasons for this. The first is that one expects finite-time singu- 
larities to form for a broad spectrum of initial data. Indeed, such singularities are 
inevitable if the scalar curvature is strictly positive. The second reason is that one 
expects the geometry of a solution to resemble a standard model (for example, a 
self-similar solution) in a space-time neighborhood of a developing singularity. The 
third reason is that having a sufficiently detailed picture of a developing singular- 
ity facilitates the geometric-topological surgeries by which Ricci flow decomposes a 
given manifold. 

Whenever a compact solution (M™, <?(•)) of Ricci flow encounters a singularity 
at time T < oo, standard short-time existence results imply that 

lim max |Rm(x,i)| = oo. 

t/°Tx£M n 

Recent results of N. Sesum (2HI allow one to replace this condition by 

lim sup max | Kc(x,t)\ = oo. 

The most interesting cases are those where a local singularity forms, that is, where 
there exists an open set fi C M n such that 

sup | Rc(a;, t)\ < oo. 

fix[0,T) 

The first rigorous constructions of local singularities were done by M. Simon 29 . 
Here the manifold is a noncompact warped product KxjS", and a supersolution 
of the Ricci flow PDE is used to prove that / vanishes in finite time on a proper 
subset of R. A second class of examples was constructed in ■ Here the manifold 
is a complex line bundle C <— ► L" fc — » CP"~ with 1 < k < n— 1. As the singularity 
forms, the flow performs an algebraic-geometric blow-down of the CP™ -1 , while the 
rest of the manifold converges locally smoothly to a metric cone on (C n \{0})/Zfe. 

As G. Perelman writes, "The most natural way of forming a singularity in finite 
time is by pinching an (almost) round cylindrical neck." [23 This is the kind of local 
singularity that we analyze in this paper. It is a continuation of our earlier work 0, 
where we gave the first rigorous examples of local singularity formation on compact 
manifolds by constructing neckpinches for rotationally-symmetric metrics on S n+1 . 
In that paper, we obtained local a priori estimates for the space-time scales at which 
a developing neckpinch singularity resembles the self-similarly-shrinking cylinder 
soliton. The present paper shows that those estimates are sharp, and gives precise 
asymptotics for neckpinch formation. 
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When considering a local singularity, two natural and important questions arise. 
(1) What is the nature of the set of points in space on which the metric becomes 
singular? (2) What is the asymptotic behavior of the solution near this set as 
the singularity time is approached? In the past two decades, a rich literature of 
both rigorous and formal matched asymptotics has developed for analyzing the 
local behavior of singular solutions of nonlinear pde such as u t = A (log u) or 
u t = Au + F(u), where F(u) = u p or F(u) — e u . Some of the many noteworthy 
results of this type are [32], HJ, QU, Q21 El C3 , H3E1E2, H3, and Q3j. 
A few results specific to geometric evolution equations such as mean curvature flow 
and harmonic map flow are PP, @][5], 0j an d 0- 

1.2. Overview. This paper is divided into two main parts. 

Rigorous asymptotics. In Section [21 we derive rigorous asymptotics for rota- 
tionally symmetric neckpinches. The set-up is as follows. Let g caxl denote the round 
metric of radius 1 on S n . Then any SO(n + l)-invariant metric on S n+ can be 
written as 

g = ip 2 dx 2 + ijj 2 g cml 

on (—1,1) x S n , which may be naturally identified with the sphere S n+1 with its 
north and south poles removed. The quantity ij){x, t) > may thus be regarded as 
the radius of the hypersurface {x} x S n at time t. It is natural to write geometric 
quantities related to g in terms of the distance s (x) = J* ip(x) dx from the equator. 
Then one can write the metric in the nicer form of a warped product 

g = ds 2 + -0 2 #can, 

bearing in mind that s is ultimately a function of both x and t. (Note that one 
encounters the commutator [dt, d s ] = —nip~ 1 ip ss d s when taking derivatives with x 
held constant.) 

In 0, we established neckpinching for a class of data essentially described by 
three conditions: (1) the initial metric should have positive scalar curvature on 
all of S n+1 and positive Ricci curvature on the polar caps (in the terminology 
of P]); (2) its sectional curvature should be positive on planes tangential to the 
spheres {x} x S n ; and (3) it should be 'sufficiently pinched', i.e. the minimum 
radius should be sufficiently small relative to the maximum radius. In this paper, 
we impose an additional hypothesis: (4) the initial metric should be reflection 
symmetric, i.e. tp(s,0) = ip(—s,0). Under these hypotheses, the results we obtain 
in Section [21 (combined with results from the predecessor [3] to this paper) imply 
the the solution obeys a precise asymptotic profile. 

To describe this profile requires some additional notation. If T < oo is the 
singularity time, let u = ip / \j2(n — 1)(T — t) be the blown- up radius defined below 
in (J2J) , let a = s / \JT — t be the rescaled distance to the neck defined in ; and let 
t = — log(T — t) denote the rescaled time variable defined in JSJ). Note that (er, r) 
are self-similar coordinates with respect to the blown-up flow. We summarize our 
results as 

Theorem 1. For an open set of initial metrics symmetric with respect to reflection 
and rotation, the solution (S n+1 , g{t)) of Ricci flow becomes singular at T < oo. 
Its diameter remains bounded for all t G [0,T). The singularity occurs only on the 
hypersurface {0} x S n . The solution satisfies the following asymptotic profile. 
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Inner region: on any interval \a\ < A, one has 

a 2 -2 1 

u(cr, t) = 1 H h o( — ) uniformly as t — ► oo. 

8t t 

Intermediate region: on any interval A < \o~\ < B^/t, one has 



u(o~, t) = yl + (1 + o(l)) — uniformly as r — > oo. 
V 4r 

Outer region: for any e > 0, i/iere exist C < oo and i < T such that 
Vra - 1 - e J = < -0(5,0 < f -Vn^l + e 



/or |er| > CS/r and i G (i,T). 

Formal asymptotics. In Section |3 we derive formal matched asymptotics for 
fully general neckpinches, without any symmetry assumptions whatsoever. These 
are analogous to the formal asymptotics for the m = 2 (presumably typical) case 
of mcf singularities considered in Our method here is to study arbitrary per- 
turbations of the self-similar cylinder soliton. That is, we consider g = g + h, 
where 

g = dx 2 + 2(n - 1)(T - t) g can 

is the cylinder soliton on R x S n , and h is an arbitrary (small) (2, 0)-tensor. To 
accomplish our analysis, we modify the Ricci flow so that the soliton becomes a 
fixed point of a related strictly parabolic flow. Because the linearization of that flow 
has a null eigenvalue, one expects to find a center manifold. We proceed by carrying 
out a quadratic variational analysis in order formally to compute the dynamics of 
the flow on that center manifold. 

The computations involved are extensive, but the conclusion we obtain is quite 
satisfying: the formal matched asymptotics suggest that the behavior analyzed rig- 
orously in Section[5]should indeed be stable for fully general Ricci flow neckpinches. 



1.3. A diameter bound. We conclude this introduction with an observation that 
follows directly from the estimates obtained in This result confirms an expec- 
tation of G. Perelman |19| . 

Lemma 2. Let (S n+1 , g(t)) be any SO(n + \)-invariant solution of Ricci flow such 
that g(0) has positive scalar curvature and positive sectional curvature on planes 
tangential to the spheres {x} x S n . 

Assume the metric g(t) has at least two bumps for all t < T (in the language of 
Let x — a(t) and x — b(t) be the locations of the left- and right-most bumps, 
and assume that for all t < T one has tp(a(t),t) > c , ip{b{t),t) > c for some 
constant c > 0. 

If g(t) becomes singular at T < oo, then diam(S' n+1 , g(t)) remains bounded as 
t/T. 
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Proof. By Proposition 5.4 of 0, the limit profile ip(-,T) exists. Let a(t) < b(t) in 
(—1,1) denote the left-most and right-most bumps, respectively. Then a(t) —> a(T) 
and b(t) — > b(T). By Lemma 5.6, the Ricci curvature is positive (so distances are 
decreasing) on (—1, a(t)] and [b(t), 1). Hence it will suffice to bound dist g r t \(xi, x 2 ) 
for arbitrary x\ < x 2 in (o(T) - e, b(T) + e) C (— 1, 1). 



Equations (5) and (11) of 3 imply that 

f-r-2 




Proposition 5.1 of [H] bounds tp a uniformly, while Lemma 5.5 shows that the number 
of bumps and necks is nonincreasing in time. It follows that 

"" (ti\ 2 ds < c r M < cl ' - » ] < T ^ 7T . 

sOi) V V / isfxi) V V'romC*) Vmax\t) VroinW 

Hence Lemma 6.1 lets us conclude that 
d 



•dist ff ( t) (xi,x 2 ; 



c* c 
< — < 



which is integrable. □ 

Recalling Lemmata 7.1 and 10.1 of [3], one immediately obtains 

Proposition 3. Assume that (S n+1 , g(t)) satisfies the hypotheses of Lemma\^and 

is reflection-symmetric. Then the final-time profile ip(-,T) has ij){x,T) > for all 
i^O. Because \ Rm | < C/ip 2 , this implies that the singularity occurs only on the 
hypersurface {0} x S n . 



Diameter estimates for arbitrary finite-time Ricci flow singularities are studied 
from another point of view in |3(J) . Lower bounds for diameter are known for certain 
topologies; see 

2. Rigorous asymptotics of the SO(M" +1 ) invariant neckpinch 
2.1. The blown-up radius and the linearized equation. The radius tp satisfies 

I -ib 2 

(1) & = V.. -(»-!)— 7 —■ 

We will assume that the solution is defined on some time interval to < t < T, where 
the initial time to ncccd not be 0. In order to construct examples of solutions that 
follow our precise asymptotics, we will at some point actually have to assume that 
T — to is "sufficiently small" . 
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We consider the blown-up radius 



(2) u -lT77v where Rn{t) = ^/2{n-l){T-i), 

H n (t) 

in which T may or may not be the blow-up time. We also introduce the rescaled 
distance to the neck 



(3) 



(Here the absence of a factor n — 1 in \/T — t is intentional and will result in the 
numerically simplest equations later on.) Then we have 



ips = y / 2(n-l) u a , %l) ss = y/%p~V)—j==, 

and 



= y/2(n - 1)(T -t)u t - ^^—^ u. 

This leads to the following evolution equation for u, 

1 / 1 \ u 2 

(4) Ur = u arT + - «-- + ( n -l)-£ 

2 \ u J u 

where 

(5) r = -log(T-t). 

The blown-up time variable r takes values in the interval 

r < r < oo, where To = — log(T — to). 

Without loss of generality we will assume that r > 1. 

Near the neck, u will be close to the self-similar radius u s=s 1. So to linearize, 
we set 

(6) u = 1 + v 
and observe that v must satisfy 

(7) W T =tW+V + - -i-2 1_. 

1 + V 

Here the partial derivatives d/da and <9/9t do not commute. One has instead 

[d t ,d s ] = -n^-ds, 

and thus 

(8) [d T ,d a }= [(T-t)d t ,Vr~td s ] =-|i + n ^|a <r . 
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2.2. Pointwise estimates for u and v. In we showed that on solutions of 
the Ricci flow whose scalar curvature is nonnegative, the quantity 

F=^2 + log. L 



L { °L min (0) 
satisfies a maximum principle. Here 

K = — and L = — — . 

and 



imin(O) = inf 



t=0 



ip 2 

The maximum principle implies that sup F does not increase when it is above n — 1. 

Let us assume throughout this paper that the initial metric has nonnegative 
scalar curvature, and that 

supF < F* 
for some fixed constant F* ^ n — 1. 

When one writes the quantity F in terms of ip, ip a , and ip ss , the estimate F < F* 
implies the following upper bounds for ip, tp a , and ?pss in terms of the radius of the 
neck i/)(0,t). 

Lemma 4. There exist constants S > and C < oo that only depend on F* and 
the dimension n such that 

C 1 + rr 2 

(9) 1 < u < 1 + C— — for \a\ < 2^r, 

r r 

and 



(10) 1 - - < u < /log M /or 2 ^r < | CT | < e s \ 

T v T V V T 

For tte derivative u a , we have 

(11) Kl < C*i^M / or | CT | < 2% /r, 



(12) | Uff |<_^ : ^/i og M /or 2Vr<| ( 7|< e ^. 

In both cases, these estimates imply that \u a \ < CS at all points with \o~\ < e Sr . 
After decreasing 5, if necessary, we may assume that CS is as small as we like. In 
particular, in the region we are considering, we always may assume that 

(13) K| < \- 

Since v = u — 1, we immediately get estimates for v and v„. Those for v a are of 
course the same as for u a . For v, we have 



(14) 



- < v < C^-^- for H < 2^7, 

T T 



-<v< cMJlogM for 2^r < \a\ < e*\ 
T V r V V T 
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in which C is again a constant depending only on F* and n. 

Proof of Lemma^ The upper estimates follow immediately from Lemma 9.4 in 
Only the lower estimate in (|5J| was not proved in pj]. Briefly, the lower estimate 
follows from the boundedness of F. At the neck, one has 

til L - — 

ip ' ip 

so that {K/L) logL < C implies 



K = Z£, L=—, ^ = (l + o(l))V2(n-l)(r-t), 



< Ipss < — 7- 
Tip 

The radius ^(0, t) of the neck satisfies 

71-1 n-l-Cir- 1 ) 
ip ip 

which upon integration yields 

mt) > ^1 - ^ y/2(n-l)(T-t) 

where C only depends on F* and n. This implies the lower estimates in the Lemma. 

□ 



2.3. Transition to commuting variables. On most of the neck, the quantity 
v = oil). So it is natural to drop the quadratic terms in Q, resulting in the 
deceptively simple equation 

Vt = Vera + V 

for v. Here the partial derivatives d T and d a do not commute. To work with 
commuting variables, we regard v as function of a. In other words, we define 

1 



(15) 



a(x, t) 



ds 



-r/2 



ds , 



Vr~tJ Q 

and consider a function v{o~,t) such that v(o~(x, r), r) = v(x, r). We shall abuse 
notation and write v for both quantities v and v. This is only ambiguous when we 
take time derivatives, in which case we write 



dv 

The two are related by 



dv 



dv dv 
= and ■F" 

OT OT 



dv 
Or 



dv 
dr' 



da dv 
dr da 



Lemma 5. Assuming reflection symmetry of the metric, i.e. ip(—s,t) = ip(s,t), 
one has 

da a . . 
- = -+nJ(a,r), 

where 

,,2 



J{a,r) 



o 1 + w 



do - 



i+f y (i+w) 2 



da. 
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Proof. By definition one has a(x, r) = e T ^ 2 s(x, T — e T ). Hence 

a T = icr + e T/2 St{x,T - e" T )e" T = ^a + e" r/2 s t (x, T - e~ T ). 
One also has (31 §10] 

8-t=J — ds = n {- + J ^ d3 \- 

Using the relations (0) and © between u = 1 + v and and s and a, one finds 
the stated expression for d T a. □ 



From here on we will consider v as a function of a and r. All r derivatives are 
intended to be time derivatives with a kept constant. It follows from J7J that v 
satisfies 



cr „ . . (n-l)«J-ii 



(16) w r = v ff(T - -tv + v - nJ(a,T)v a + 

2 1 + u 

We write this equation as 

v T = A(v)+N(v), 
where A is the linear differential operator 

. f d\ 2 ad 
A= {da) -2^ + 1 ' 

and N(v) represents the remaining (nonlinear) terms in l|16fl . 

N(v) = —nv a / da J — 



o 1 + v 1 + v 



Integrate by parts to get 



(17) N(v) = -nvJ -!^da-^—±-. 

j {l + vy i + v 

Lemma 6. For all \a\ < e Sr , one has 

\N(v)\<C±±^. 



Proof. This follows from Lemma 0J In fact, Lemma 0] implies that 
(18) \v\ < C^-^- and \v„\ < 



for | a | < 1\Jt. For larger a, one finds other estimates for v and v a , which are 
stronger than (fTH|) . Using (|T%)l for all \a\ < e 5r , one arrives at our estimate for 
N(v). □ 
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2.4. The linearized equation. If we ignore the nonlinear terms in (|16|l . then we 
see that the small quantity v satisfies v T — A(v). Even though v is not defined for 
all s S K, we know from our estimate on the neck that v is defined for |er| < e Sr , 
for some S = 5{F* , n). In the first analysis, we assume that v is in fact a solution of 
the Cauchy problem v T = A(v). The operator A is self-adjoint in the Hilbert space 



H = 



(v € L 2 (m^e-^^da^ \ v(a) = v(-a)\ 



It has pure point spectrum, with eigenvalues A m = 1 — m for m = 0, 1, 2, . . .. The 
eigenfunction corresponding to the eigenvalue X m is an even Hermite polynomial 
of degree 2m. Normalizing the eigenfunction so that the coefficient of its highest 
order term is 1, we set 

h (rA ^2m ( 2m V- 2m-2 , ( 2to ) ! 2m-4 

/l2m(CT)=a -l!(2m-2)! CT + 2!(2m-4)!" _ '" 
In particular, 

(19) h {a) = l, h 2 {a) = a 2 - 2, h 4 (cr) = ct 4 - 12a 2 + 12. 



2.5. Eigenfunction decomposition of u. Let r\ E C°°(R) be an even bump 
function with i](z) — for z > 1 j and = 1 for z < 1. We then define 

(20) F(cr, r) = ?7(CTe"' 5r/2 )w(CT, r). 

For |er| > fe 5r / 2 , we set V = 0. Then V(-,t) € W, and w = V for |cr| < e Sr ' 2 . 
A computation shows that 

(21) V T - A(V) = r]N(v) + E, 
where the "error term" is 

E = (rj T - T\aa + T^Virj v - ^VcrVa- 

Lemma 7. The error term E vanishes except when e Sr / 2 < \a\ < je Sr ^ 2 . When 
E ^ 0, one ftas 

|£(<r,r)| <C|(t|, 
where the constant C only depends on F* and n. One also has 

\\E(;T)\\<Cexp(-e ST / 2 ), 

where the constant C again only depends on F* and n. 

Proof. The pointwise estimate follows from the boundedness of r] T , 0~r\ o , and r] aa , 
as well as the pointwise bounds for v provided by Lemma 01 Given the pointwise 
bounds for E, one finds that 

/>oo 

\\E\\ 2 <C 2 o- 2 e-° 2 / 4 da, 

J e Sr/2 



from which the TL norm estimate follows. 



□ 



ASYMPTOTICS OF THE RICCT FLOW NECKPINCH 



11 



For each r, the function V(-,t) belongs to TL. So we can consider the splitting 
into mutually orthogonal terms determined by 

(22) V{a, T ) = a (T)h (a) + a 2 (T)h 2 (a) + W(a,r). 

Our pointwise estimate 114JI for v suggests that for large r these terms will decay 
like t -1 . This turns out to be the case for the middle term a2(r)h2(r), but the 
other two terms are in fact smaller. We will show in the next few paragraphs that 
they decay like r~ 2 . This establishes 02(1") /i 2 (t) as the dominant term for large r. 



2.6. Easy estimates for ao and a 2 . Lemma 1 and in particular (|14|) imply that 



\V(a,r)\<-(l + a') 

for all cr G K, and all t > tq. Upon taking the 7i-inner product with ho and h 2 this 
gives us 

(23) |oo(r)| + Mr)|<! 



for all t > tq. 



2.7. Decay of a . Since ||/io|| 2 ao(r) = {h$,V)u, we have 

\\h fa' (T) = {h ,V T ) n 

= (ho,A(V)+ V N(v)+E)n 
= {h ,V)n + {ho,vN{v)+E) n 
= \\h \\ 2 a (r) + (ho,vN(v)+E) n . 

Hence 



a' (r)-a (r) = (j^,r,N(v)+E\ 

Using our estimates for N(v) and E, we conclude that 

fo(r) #ao( T ) -«o(r) 

satisfies 

IAMI < C(F *'" ) 



n 



T 

The variation of constants formula tells us that 



oo(r)=-y e T - r / (r')dr', 

and thus that 

. C(F*,n) 
\a (r)\ < V ' 
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2.8. Decay of W. Since W lies in the stable space of the operator A, i.e. the space 
spanned by those eigenfunctions with negative eigenvalues, we expect W to decay 
according to the slowest stable eigenvalue (i.e. ~ e _T ), or else to decay like the 
"forcing term" rjN(v) + E. In this situation, the forcing term is dominant, and we 
have the following. 

Lemma 8. For any given tq > 0, one has for all r > tq 

\\W(;T)\\<e^\\W(;To)\\+^. 
Consequently, one also has 

\\W{;r)\\<^(l + W ). 
Here the constant C only depends on F* and n, and 

W =e To \\W(;T )\\. 

Proof. Since \\W\\ 2 = (W,W) = (W, V) and {W,A(V)) n = {W,A(W)) H , one has 

= (W,V T ) n 

= (W,A(V) + nN(v) + E) n 

= (W, A(W)) H + {W, V N(v) + E) H 

< -\\W\\ 2 + \\W\\ ■ \\j]N{v) +E\\, 

whence 

(24) ^D<_||w|| + ||^(t;)+f?||. 

dr 

From Lemma|Hlwe conclude that \t]N(v)\ < C(l + er 4 )r~ 2 for all cr, and hence that 

II^WII < 

For E, we use Lemma together with the calculus inequality Ce~ e ' < C't~ 2 . 
We apply all this to (|24|l and conclude that 

Therefore 



d{e T \\W\\) 
dr 

Integration leads to 



< Ce 



T T- 2 . 



e T ||^(-,r)|| < e ro ||^(-,r )||+C / Zrdr <e T °\\W(;T )\\ + S-e\ 



T 2 

TO ' 



as claimed. The second estimate for ||W|| follows directly from supr 2 e T < oo. □ 
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2.9. Derivative estimates for W. We use the regularizing effect of the heat 
equation to bootstrap the estimates of Lemma |H1 by one space derivative. 

Lemma 9. Let r > be as in Lemma\^ Then for all r > r + 1, one has 

||W CT (-,r)|| <C(1 + W )t- 2 , 

where the constant C only depends on F* andn, and as before, Wq = e T ° \\W(-, tq)\\ . 



Proof. Note that 



Au — u aa — 2 (JU(7 + u — e a ^{e a / 4 it 



so that for u G D{A) one has 

-(u,Au) H = [ e-° 2 / 4 {ul-u 2 }da. 
Jn 

Thus 

(25) HM| a = \\u\\ 2 -(u,Au). 
If u _L ho, then (u, — Au)-h > 0, and we have 

\\dau\l < \\u\\ + y/(u, -Au)n- 

Applying the spectral theorem to the selfadjoint operator A on the Hilbert space 
{ho, /12} 1 " in order to estimate (u, Ae eA u), one concludes that 

(26) \\dy A u\\ < ^-\\u\\ 

holds for < 6 < 1, and all u S D(A). 

The stable component W of V satisfies 

W r = A{W) + /(r), where /(t) = T(nN(v) + E), 

and where V is the 7i-orthogonal projection onto the space {u G Tt \ u _L ho, h^}- 
By the variation of constants formula one then has, for all r > tq + 1 



W{t) = e A W{T - 1) + / e 0A f(r - 8) d6. 
Jo 

Since f(r) _L {ho, /12}, we can use to get 

RW(r)|| < \\d a e A W{T-l)\\+ C \\d a e eA f{T-6)\\ d9 

Jo 

< C\\W(t - 1)|| +J °\\f(r- 9)\\ d6. 

The combination of Lemmas El and [7| again gives ||/(t)|| < Ct~ 2 . The integral 
can therefore be bounded by CV -2 . Combining this with our previous estimate for 
||H^(t — 1)|| from Lemma 03 we get the inequality in Lemma [5] □ 
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2.10. A bound for ||oW||. If <p € C^(R), then one has the identity 

^ a2 / 4 2<p(a)<7cp / (a)d<j= f e- a2 / 4 (-2ip(a) 2 + ±a 2 p{a) 2 ) da 
Jn 

which after rearranging leads to 

\\<J<p\\ 2 = 4 ( e- a2/4 (a<p<p' + p 2 )da 



= 4{(a<p,<p') H + \\<p\\ 2 } 

<i||H| 2 + 8||^'|| 2 + 4MI 



<i6{y\\ 2 + y\\ 2 }. 

An approximation argument leads to 

Lemma 10. For any ip € TL with ip a e TL one has op e "H, while 

\\a V {a)\\ <4{||^|| + \\<p\\} . 

In particular, W satisfies 

(27) \\aW\\ < C(l + W )t- 2 

for all t > To + 1 . 

2.11. Pointwise estimates for W . A Sobolev inequality says that the bound on 
1 1 o'er W 1 1 implies a pointwise bound for W. Due to the exponential weight in the 
norm this bound is not uniform. One has 

Lemma 11. Any ip £ H with d a p e TL is a continuous function, which satisfies 

e-° 2 / 8 \ v (<j)\<c(M + \\dM)- 

Hence W satisfies 

m *,r)\< C ( 1+ 2 W °h^. 

Proof. Assuming p 6 C*(R) one has 

e V/ V(a) 2 = ^ e^ 2 ^[~^ 2 + 2 Wa ]da 
J —oo 2 

< \\<r<p\\\\<p\\ + imiikii 

<C{\\p\\ 2 + \\pJ 2 }. 
So the pointwise estimate holds for p £ C 1 . By approximation it also holds for all 

(pen with d a tp en. □ 

Lemma 1771 savs that W decays pointwise like t -2 , and uniformly on bounded a 
intervals. The following pointwise bound for W only gives a r _1 decay rate, but it 
is stronger for large values of a. The estimate also does not depend on Wq. 

Lemma 12. For some constant C = C(F* , n) < oo one has for all r > tq + 1 

\W(a,r)\<-(l + a 2 ). 

T 
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Proof. We have, by definition, W(a, r) = V(a, r) — ao(T)ho(a) ~ a2{r)h2{o-). Using 
our pointwise bounds for V, ao, 02, and ho(a) = 1, /12(c) = a 2 — 2, one quickly 
gets the stated estimate for W. □ 

2.12. Pointwise estimate for W a . By definition we have W(t, a) = V(t, a) — 
ao(r)ha(o~)—a 2 {T)h 2 (o~). Differentiate and keep in mind that h' (a) — 0a,ndh 2 (a) = 
2a. We get W a — V a — 2a2{r)a. Using the pointwise derivative bounds and 
(|T2"j) for u, one then gets 

(28) |wg<£(i + H ) 

which holds for all t > tq and a £ R. 

2.13. An equation for 02- As with ao, we have 

\\h 2 \\ 2 a' 2 {r) = {h2,V T ) n 

= {h 2 ,A(v) + r]N(v) + E) n 

= {A(h 2 ),V) H + (h 2 ,r]N(v)+E) H 
= {h2,nN{v)+E) H . 

We use (|17l) and the identity 1/(1 + v) = 1 — v/(l + v) to rewrite N(v) as 
(29) 

N 2 {v)=-{vl + \v*) 

N(v) = N 2 (v)+N 3 (v), where { r° J wl+hv 3 

N 3 {v) = -nv a - ° da+ ! - 



{1 + v) 2 1 + v 

Here, N 2 (v) is the purely quadratic part, while N 3 (v) contains the cubic and higher 
order terms. 

Lemma 13. There is a constant C — C{F*,ri) such that for all r > tq 

WvN 3 (v)\\ < ^. 

Proof. The pointwise estimates for v and v a from 1|14|) imply that 

T T 



whenever \a\ < e Sr . Hence we have 

1 1 _6 



.V,(c) _ C^-^ for \a\ <_.<" 



T 

(The term with v 3 contributes the highest power in a.) 

This implies the Lemma. □ 

We continue with our computation of a' 2 (t) . We have 

\\h 2 \\ 2 a' 2 (T) = (h 2 ,r 1 N(v)+E)- H 

= (h 2 ,vN2(v)+vN3^) + E) H 

= {h 2 ,r)N 2 (v)) H + {h 2 ,r)N 3 (v) + E) H . 
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The last term satisfies 

\(h 2 ,r)N 3 (v) + E) H \< \\h 2 \\ (||t?JV 3 («)|| + \\E\\) < 

We would now like to replace the v in the quadratic expression N 2 (v) by V, and 
then by the dominant term in l|22[) . i.e. a 2 (r)h 2 (a). Before we do this, we estimate 
the errors produced by these replacements. 

Lemma 14. There is a constant C — C{F* , n) such that for all r > r 0; one has 

\\v(N 2 (v)-N 2 (V))\\<Ce- e6T/2 . 

In other words r){N 2 (v) — N 2 (V)) satisfies the same super-exponentially small 
estimate as the error term E. (See Lemma 0) 

Proof. Since V = r\ ■ v, a direct calculation gives 

V (N 2 (v) - N 2 (V)) - n {{ V v) a 2 - v a 2 } + ^{if - l)v 2 . 

Using the pointwise estimates for v and v a with the fact that rj (N 2 (v) — N 2 (V)) is 
supported in the region e Sr / 2 < \a\ < je Sr ^ 2 , one arrives at the stated estimate for 
the U norm of r\ (N 2 (v) - N 2 (V)). □ 

We therefore have 

\\h 2 \\ 2 H a' 2 (r) = (h 2 , V N 2 (V)) n + (h 2 , V (N 2 (v) - N 2 {V))) n + (h 2 , V N 3 (v) + E) n . 

If we write N 2 [v,w] — —v a w a — \vw for the natural symmetric bilinear expression 
with N 2 (v) — N 2 [v,v], then the first term above can be written as 

{h 2 ,vN 2 (V)) n = (Tih 2 ,N 2 (a 2 h 2 + a h + W)) H 

= (a 2 ) 2 ( V h 2 ,N 2 (h 2 )) n 

+ 2a 2 {rjh 2) N 2 [h 2 ,a h + W]) H 

+ (vh 2 ,N 2 (a h + W)) n , 

which is a quadratic polynomial in a 2 . Adding in the omitted error terms, we reach 
the following observation. 

Lemma 15. 

4(r) = K(T)a 2 (T) 2 + 2L(r)a 2 (T) + M(r), 

in which 

K(t) = jXr 2 {vh2,N 2 {h 2 )) H , 

L ( T ) = T1A2 (vh2,N 2 [h 2 ,a h + W)) H , 

M(r) = {(r,h 2 ,N 2 (a Q h Q + W)) H + (h 2l r, (N 2 (v) - N 2 {V))) H + {h 2 , V N s {v) + E) H } . 

Lemma 16. 



K(r)=-8 + 0(e~ e6T/2 ) 
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Proof. Since N 2 (v) = — v 2 — ^v 2 , the explicit expressions i(T5|) for h , h 2 , and /14 
give us 

N 2 (h 2 ) = -^cr 4 - 2(T 2 - 2 = -|/i 4 -8/12 - 12 /i . 
Hence (h 2 ,N 2 {h 2 ))n = ~8\\h 2 \\ 2 . This implies 

^ \ (h 2 ,N 2 (h 2 )) n ^ ((l-y)h 2 ,N 2 (h 2 )) n 

K{T) = w 5 + m a ■ 

The first term is —8. The other term can be written as the integral of a function 
with fixed polynomial growth, which is supported in the region \o~ I > e ST ' 2 . This 
leads quickly to the stated estimate. □ 

Lemma 17. L(t) = O (t~ 2 ), i.e. there is a constant C — C(F*,n) such that for 
all t > To + 1 one has \t 2 L(t)\ < C . 

Proof. Using (bi)lincarity of N 2 [v,v], we get 

\\h 2 \\ 2 L{r) = {nh 2 ,N 2 [h 2 ,a h + W]) n 

= ao(r) (h 2 ,riN 2 (h 2) ho)) n + {h 2 ,r)N 2 (h 2 ,W)) n . 
The first term is C(r~ 2 ). For the other term, we have 

\(h 2 ,r)N 2 (h 2 ,W)) n \ = \ {vh 2 ,ti 2 {a)W a + ±h 2 W) H \ 
< \\vh 2 h' 2 \\\\W a \\+C\\W\\ 
C 

by our estimates for \\W\\ and ||9 CT W||. □ 
Lemma 18. 

M(t)<C\ || W (. )T0 )|| + - 

^ T 7 

In fact one can show that M(r) < C Q T~( 4 ~ a ' for any a > 0, but we will only 
need the estimate with a — I. 

Proof. In Lemma El ||/i2|| 2 M(r) is defined as the sum of three terms. We now 
estimate them one by one. 

The middle term is 

\(h 2 ,ri(N 2 (v) - N 2 (V))) H \ < \\h 2 \\\\r](N 2 {v) - N 2 (V))\\ <Ce- &&Tl \ 

by Lemma ITU 

The third term is bounded by 

\(h 2 , V N 3 (v) + EU<C(\\nN 3 (v)\\ + \\E\\)<^, 

because of Lemma |l 31 and the super-exponential estimate for E in Lemma 

The first term can itself be split into three terms: 
(30) 

( V h 2 , N 2 (a h 2 +W)) H = ao(r) 2 (^2, N 2 {h 2 )) H +2a (T)( V h 2 , N 2 [h 2 , W]) H +(r]h 2 ,N 2 (W)) H . 
The first of these terms is 0(r~ 4 ), since ao(i") = 0(t~ 2 ). 
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From 



{vh 2 ,N 2 [h 2 ,W]) H = (r)h 2 , -h' 2 (a)W a - \h 2 W) n 

1 

WW) 



= - ( V h 2 h' 2 , W a ) H -^- ( V h 2 (a) 2 , W) H 

= 0(WWJ 



and ao(r) = 0(t~ 2 ), we deduce that the middle term in 1)30(1 is also 0(t~ 

We are left with the last term in <|3(J|I . This term can be expanded as 
(31) (r)h 2 ,N 2 (W)) n = -(nh 2 ,W*) H - (i!h 2 ,\W 2 ) H . 

The second term satishes 

?- a2 ^rih2(v)W(T,a) 2 d<T < C(||M^|| 2 + \\<rW\\ 2 ) 

<c{e 2 ^\\W(;T Q )\\ 2 + ± 

by Lemmas 151 and ITU1 

We bound the first term in (|31fl as follows: 

\(vh 2 ,W ( 



r <?)n\ 



e " ' 'r)h 2 (a)W a da 
f e-° 2 l\l + a 2 )W 2 a d<j 



< c 

< CWil + afWJlWal 



The pointwise estimate (|28|l bounds the first factor by 



||(l + a) 2 W CT || < 



G 



while Lemma implies 



\\W„\\<C e ro - T ||W(-,r )|| + 



Combined, these last two estimates imply Lemma IT5I 

2.14. Dichotomy for the decay of ci 2 (t). We have shown that 

da 2 (r) 



(32) 
where 

Hence 



dr 



-8a 2 (rY+M*( T ), 



\M*(t)\ < C(1 + W )t- 



r d{Ta2 j T)) = (ra 2 (r)) - 8(™ 2 (r)) 2 + r 2 M*(r). 
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Lemma 19. The coefficient 1x2 (t) satifies either 

(33) 02 (t) <C(1 + W )t- 2 
/or oZZ t > ro + 1, or e/se 

(34) lim ™ 2 (r) = -. 

r — >oo o 

2.15. Pinching time. We will show below that alternative (|34f) does indeed occur, 
namely that there exists an open set of solutions matching the asymptotic profile 
we have constructed. Our first step is to improve the estimate obtained in Lemma 
6.1 of 0. Let r(t) = t/j(0, t) denote the radius of the neck, and recall that 

K=^f and L = Llf 

denote the sectional curvatures. Define A = 2 — logL min (0), so that we may write 
the cylindricality F as 

F = £ (A + logi) = _^( A + log I_*l 

Recall that there exists F* > n — 1 depending only on g(0) such that < F < F* 
on the neck. Writing the evolution equation for the radius 



ip ss - (n- 1)- 



= Ips 
s=0 



dr 
~dt 

in the form 

(35) -§(r 2 ) = -2(»-l) + 2 f i 
at A + log ^ 

then proves that 

(36) -2(n- 1) < ^-(r 2 ) < -2(n- 1) + 2 



n- 1 



2 ) ^ -_:!/ ■ - 1. I T ■ p. 

A + log 



This observation lets us prove that the singularity time depends continuously on 
the initial conditions. 

Lemma 20. For rotationally symmetric metrics g(t) that develop a neckpinch, the 
pinching time T is a continuous function of the initial metric go ■ 



Proof. Given a metric go, let To = T(go) denote the singularity time of the solution 
g{t) with initial data go- 
To prove lower semicontinuity, let T_ < To be given. Then {g(t) : < t < T_} 
is a smooth solution of Ricci flow. Because regular solutions of parabolic equations 
depend continuously on their initial data, a solution g(t) will remain smooth for 
<t < T_, provided that its initial data go is sufficiently close to go- This implies 
that the singularity time T of g(-) satisfies T > T_. 

To prove upper semicontinuity, let T + > To be given. At any time T* < To, 
observe that for all go sufficiently close to go, the radius f of the solution g(t) 
satisfies f(T,) 2 < 2r(T*) 2 , and its cylindricality F satisfies F < F* < 2F* . By 
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taking T* close enough to To, one may by l|36(l ensure that any solution g(t) with 
initial data go satisfies 

d 4F* 
T Ar 2 ) < -2(n - 1) + - < -(n - 1) 

at A + log 2r ^yi 

for all t € (T* , T) , where f is the pinching time of g(t). So for t g (T* , T) , it follows 
that 

f(t) 2 <f(T*) 2 - (n-l)(t-T*) 

< 2r(T,) 2 - (n-l)(t-T*) 

<4(n-l)(T -T»)-(n-l)(t-r»). 

Here we used the fact that r(t) 2 < 2(n - l)(2b - Choose T* < T closer to T 
if necessary so that T* + 4(T — T*) < T + . Then if t > T* + 4(T - T*), one has 
f (t) 2 < 0. This implies that f <T. Jf + 4(T - T*) < T+. □ 

2.16. Initial data. Given a target time T, we now construct initial data u = 
1 + u at t = to <T that will become singular at T and satisfy the asymptotic 
profile (|34() . In the next subsection, we shall obtain appropriate estimates for the 
function V(ct,t) = -q{e~ ST / 2 o) ■ v(ct,t) = X)^o fl2fe ( r )' l2fc ( cr ) corresponding to this 
construction. 

Let £ g (0, ^-T=) be a constant to be chosen later. Near the developing neck- 
pinch, we impose the initial profile 

U = 1 + a h (a) + a 2 h 2 (a), 

where «o and a 2 > are constants to be chosen later, and cr(s,t) — sj \JT — t. 
That is, we prescribe tp(s,to) = R n {to) ■ u(a,to) as 

,2 



S 

(1 + a - 2a 2 ) + a 2 -, 



(37) V(Mo) = V2(n-l)Vr-i .„ „ r , 

on the region -A/ = {s:|s|<so}, where So is given by 

1 



so 



eVr-tolog 



T-to 



We will require 

(38) -e < a Q < 2a 2 
and 

(39) — <a 2 < — , 

To t 

where 

T0=1Og T 1 ^ >>L 
The radius of a neck of the form (|37(l is 



(40) r = V(0, to) = ^2(n - 1)a/T - t (l + a - 2a 2 ). 

Thus 



a 2 s 2 



^ = r + V2(n^l)-^__. 
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Since 



we have |^ s | < ip s (so,to) onAf, where ip s (so,to) = 2ey/2(n — 1) • ro«2 satisfies 

£ 2 < ^(ao.to) < g 
" 2V2(n-l) " 

Outside A/", we extend V 1 linearly with slope ip s (so, to) > until we reach the height 
■0 = 1, where we glue on a standard profile that is smooth at the poles. As in §8 
of [3], this can be done with \i/j s \ < 1 everywhere. Notice that the initial profile 
ip(-,t ) is smooth except on dAf; the corresponding solution of the Ricci flow will 
be smooth for all t > to that it exists. 

We next obtain a formula for the upper bound F* on the cylindricality F of a 
neck of the form l|37l) . where ao and «2 satisfy (|38l) and l|39[l . respectively. Note 
that F < off Af. On Af, we have 



^ < ^(*o, h) = r + y/2(n-l)- 

V 1 ~ *o 

< V2(n-l)(l + e 2 r )^T-t a 

and 

/ o /of ~ "2 , 2^/2(n- 1) 

SS = 2V2(n - 1)-=== < /?7 ^— ■ 

Independent of e < j, we may choose close enough to T so that for any to £ 
[i» , T) , the radius is bounded from below by 



2 \ ^_ y/2(n-l) 



r > y/2(n-l)y/T-to f 1 - e - — J > L y/T-t , 

and so that the sectional curvature L is bounded from below on Af by 

r _ l-ips 1 - 8 ( n - !) £2 1 2 

V(so,io) 2 " 4(n-l)(r-to)(l + e 2 r ) 2 " e ' 

provided tq > Tq (n, e). Here we used the fact that lim x \^o(a; log -) = 0. By Lemma 
9.1 of 3 , the bound L > e 2 lets us assume A = 0. Combining the estimates above, 
we conclude that there exists C = C(n,t*) such that for any to £ [t», T), one has 

K ^ ss l-^ s 2 
F = — logL = — log 



1 — tpj ip 2 

2 1 

—^Iplpss log — 

1 - 8{n — l)e 4 ro 



< 16(n- 1) (log2 + — ) 



< Ce^ro. 

Hence we may define 

(41) F* = [n - 1) + Ce 2 r . 

This formula for i* 1 * lets us adjust our constants in order to force pinching at T . 
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Lemma 21. There exist and t* depending only on n such that for any e £ 
(0,£*], all initial times to £ [t*,T), and any a 2 satisfying 13 ty . there exists some 
ceo satisfying yStfj) such that the solution with initial profile \'8T\) becomes singular 
exactly at T . 



Proof. Proposition 5.2 of [3] implies that r is monotonically decreasing. So for all 
t > to that the solution exists, one has 



r{tf <rl< 2(n - 1)(T - i ). 



As above, we fix t* < T so that (gTJ holds for all t £ [t* , T) . Then by (JSSJ , there 
exists c > depending only on n such that 

-2(n - 1) < — (r 2 ) < -2(n - 1) + 2 -. 

at to — c 

There exists C* = C*(n, t*) such that ""^f^ T " — C* £2 - Hence we may fix 
e* = e*(C*) such that for all e S (0, e*], one has 

-2(n-l)<|(r 2 )<-2(n-l)(l-e). 

By Lemma 7.1 of the solution becomes singular at T if and only if r(T) — 0. It 
follows that 

/° - < f - to < — r^-. 

2(n-l) ~ 1-e 2(n-l) 

By ||1D). this implies that 

f1 , 9 , 2 , T-<o fl+^0^202) 2 

(l + ao-2a 2 ) < - — — < • 

1 — to L — E 

Thus if ao — 2ct2 > 0, then pinching occurs at T > T. On the other hand, if 
— 1 < ao — 2ct2 < —s, then pinching occurs at T < T . By Lemma |2*H1 pinching will 
occurs exactly at T for some ao £ [2a2 — £, 2^2]. Clearly, this ao satisfies (|38|l . □ 



2.17. The dominant asymptotic profile. For t > to, let w(<r, r) denote the 
solution of the Ricci flow corresponding to the initial data l)37|l . Then 

v{a, T ) = a ho(cr) + a 2 h 2 (cr) for \a\ < a a = <j{s ,to) = er . 

For t > To, define V(a, r) = rj(e~ ST / 2 o-)v(a : r), noting that V(-, r) £Ti admits the 
orthogonal decomposition 

(42) V(a,T) = a (T)h (a) + a 2 (T)h 2 {a) + W{a,T). 

Intuitively, one expects that ao and a 2 will be close to ao and a 2 , respectively, 
at t = To and for a short time thereafter. We now make this expectation precise. 
Define U £ H by 

U = V(-, t ) - a h - a 2 h 2 . 

At t = t , one has \V a \ = \v a \ = \ip s \/^2(n - 1) < e for \a\ < a and \V a \ < 1 
everywhere. Because |ao| ! Q ! 2 £ [0, 1], there are universal constants C (which may 
change from line to line) such that 

U(a) = for all \a\ < a 

\U{a)\ < C*(1 + ct 2 ) for all | a j > a , 
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hence such that 

\\U\\ 2 <Cl (l + a 2 ) 2 e-°^ 4 da 



l\<j\><7 

< Ce~ a ?> /5 = Ce- £2r ° 2/5 . 

Because \\W\\ 2 = (W,W) = (W,U) < \\W\\ \\U\\ at r = tq, it follows that the 
constant Wo introduced in Lemma |H1 satisfies 

(43) Wo = e T ° \\W(;tq)\\ < e r ° \\U\\ < C cx P (t - ^r 2 ). 
Similarly, one has 

a 2( T o) = 7~T~7?2 (^2, V(; T )) =a 2 + —^—2 ( h 2, U) , 
IIM IIM 

and hence 

(44) \a 2 (T )-a 2 \<Ce~ £2T ^ 10 . 

We are now ready to prove that alternative l|34ll occurs. By H32[) . one has 
^-a 2 (r) = -&a 2 (T) 2 +M*{T), 

(IT 

where \M*{t)\ < C(1 + W / )t" 3 . Bv (g^ , we may make r n so large that C(l + Wn) < 
C* . Then 



-8a 2 (r) 2 - ^ < Aa 2 (r) < -8a 2 (r) 2 + % 
By (|39H and l|44|) , we may make tq larger if necessary so that 



Co £ / x 2 

~7T-< a 2 (r ) < -. 
2t r 



a 2 (r) > a 2 {ro) - — (r - r ) > ^ ( T _ T o)- 



As long as a 2 (r) 2 > ^R-t 3 , the function a 2 (r) will be decreasing. For such times, 
one has 

d 2 Co 32 Co 64 

— o 2 (t) > -8a 2 (r ) - — > --^ - > --^ 

(IT T T Q T T Q 

provided that To is large enough, and hence 

64 , , _ e 64 
2t~o 

In particular, given any C > 0, one can make To so large that 

< j-i^ : 64 1 + Wo 

a 2 [T + 1) > J > C- — 

2t r z (r + l) z 

This justifies our assumption that a 2 was decreasing for To < r < To + 1. Moreover, 
if we choose C as in Q33JI then it violates alternative l|33|) . Hence we have proved 
the following result. 

Proposition 22. For any final time T, all sufficiently small e > 0, all initial times 
to sufficiently close to T , and any a 2 satisfying ^33D > ^ ere exists cto satisfying 
$38)) such that the solution with initial profile becomes singular at T and 
satisfies the asymptotic profile 

(45) ci 2 (t) = ho(-) as r — * oo. 

8r t 
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2.18. The intermediate region. Consider a solution that forms a singularity at 
time t = T. Define u, v, etc. as above and assume that our solution satisfies l|45|l . 

We have shown that 

o- 2 -2 1 

(46) v(*,t) = — +o(-) (r ->oo), 

or r 

uniformly on any bounded interval \o~\ < A. 
Now consider the quantity 

q = u 2 - 1. 

Lemma 23. For any constants e > and B < oo, there exist A and f such that 

2 2 

whenever A < a < By/r and times t > f. 
In other words, for 1 <C cr < B-Jt one has 

We shall prove this Lemma by computing the evolution equation for q in the 
variables t and p = cr/yV, and finding suitable sub- and supersolutions. The 
evolution equation l|49|l turns out to be a small perturbation of a first order linear 
equation, which is easily solved by following its characteristics. 



2.18.1. Equation for q. Since q a — 2uu a and q aa = 2uu aa + 2m 2 ,, one has 

n — 2 q 2 

(47) q T \ x = q aa +q+— 

2 q + l 

Here the time derivative is with constant x. Using Lemma we get 

(48) q T \ a = q aa - ^rq a + q - nJ(a,r)q a + " 



2 ' y —v. - 2 g + 
where gvl^ is the derivative of g with cr held constant and where 

r& 2 

u J u z 

We now regard q as a function of the new variable p = a/^/r. Because 
one has 



(49) q T \ p - ~q PP + -q p - 1 = ~ \ Tjlp ~ nVrJq, 



t 2'" 2 g + 1 
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2.18.2. A subsolution for j49\) . We begin with some simple estimates. Fix some 
Aq 1 and recall that to is the smallest time under consideration in this paper. 
Given B > 1, we shall choose A > Aq and t\ > tq large enough to make everything 
work. The region A/^fr < p < B 2 corresponds to A < a < B 2 yfr. By Lemma 0J 
there exist Co = Cq(B) and C\ = C\(Aq) such that for all A < a < B 2 ^fr 1 one has 



a 2 



and 

Thus we have 
and 



u < 1 + C — = 1 + C p 2 < 1 + C B 4 

T 



0<u a < Co^-^- < Ci- = Ci-^=. 

T T \J T 



< u p < dp 



0< pq p = 2pu p u<C{Ao,B) 
for A/y/r < p < B 2 . Since u > 1 — C/r for C — C(n, F*), we may assume without 
loss of generality that To ^ 1 is so large that 

-^- = ^<2q 2 p <C(A ,B) 
q + 1 i'- 1 



for A/^/t < p < B 2 and r > r . Since n^/rJq p > for all < a < e 5r ', it follows 
that 

(50) P -q p - n^Jq p + - ~ 2 %- <C(A ,B). 



We introduce the linear differential operator 

-f + P ~ 

T Ipp 2 

Then (|^|) and JSUJ) imply that <j = w 2 — 1 satisfies 

KWI < 

whenever A/y/r < p < B 2 and r > tq. 

Let Q be a smooth function with support in [0, 1] such that Q(r) < r 2 /4, and 
define 

K = max | Q"(r) | . 

0<r<l 

For ri > ro, define 

(51) go(p,r) = Pl (r) 2 Q(p/pi(r)), 

with 

Pl (r)=e(— V 2 ^. 

This function is almost in the nullspace of the linear operator £; indeed, a straight- 
forward computation shows that 

^o]=- 1 -^ = - 1 -Q"( P / Pl (r)), 

so that |£[go]| < Kq/t. 
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Now for K > Kq to be determined and any 9 <E (1/2, 1), consider 

(52) q_(p,T)=9q (p,T) + -. 

T 

One then has 

£[g_] = 0£[ 9o ] + ^[r- 1 ] < A — - - < ^— - < 0, 

t r r 

so that g_ is a subsolution for all r > 7*0. Moreover, if we choose K > Kq+C(Aq, B), 
then 

(53) C^- q] < C ^V +K °- K <0 

T 

for A/y/r < p < B 2 and r > t . 
Define the region 

(54) n= Up,r) : A <p< Pl ( r ), r x < r < coj , 

where w is determined by pi(w) = £> 2 , i.e. 

w = 7i+ log (B 4 n/v4 2 ) . 

We will choose A and T\ so that .B 4 Ti > A 2 . First, we compare the values of q and 
q_ along the parabolic boundary of Q. 

Along the left edge p — A/y/r, t > t\, it follows from (|46|l that 

(55) q(A/V^, t) = u 2 {A, t) - 1 = ^— - + o(-) > (1 - E ) A 

4r t 4t 

where £o = £o(^o, Ti). On the other hand, there exists z\ = £\{A, K 1 To) such that 

6 , K , . 6 , 
9-(p,r)< i p 2 + -<(l+ £l ) i P 2 

for p > A/y/r. So q_ < g at (A/y/r, r) provided that 

i<.<i^. 
2" 1 + ci 

Along the right edge p = pi(r), one has 

M Pi (t),t)= Pi (t) 2 Q(1) + - = -, 

r r 

because Q is supported in [0, 1]. But since q p > 0, estimate (|55f) implies that 
g(pi(r),r) - q(e/ T -^' 2 A/ y/T^r) > q{A/yf?,r) > — , 

4 T 

where Sq = eq^Aq, n). So q_ < g at (/?i(t), r) provided that 



1 - e 



The bottom of is the single point (A/y/fi, r%). So the choice 8 < (1— £o)/(l+£i) 
ensures that q- < q there as well. 

Now choose K depending on Aq, B, and Kq to satisfy K > Kq + C(Aq, B). Take 
Ti > To large enough so that £o < 1/2 and choose A > Aq depending on K so that 
A > \/8K and so that e\ > is as small as desired. Then increase r\ if necessary 
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so that t\ > A 2 /B 4 , making £o as small as desired. These choices make it possible 
to take 9 G (1/2, 1) as close to 1 as one wishes and still satisfy 6 < (1 — sq)/(1 +£i). 
Furthermore, they ensure that g_ < q along the parabolic boundary of Q. Because 
e( T ~ Tl )/ 2 A/ A /Ti < B 2 for all r < oj, it follows from J53J that q > q- throughout 
Q. Since for any e, one may take Q(r) — r 2 /4 for r G (e, 1 — e), this implies in 
particular that 

2 

q(p,r)>9^ 

for all points p such that A/^/t < p < B and times r such that 

2 log {By/n/A) <r-n< 2 [log (BVnM) + log B]. 

Our estimates only improve if we increase t\ (without altering any other choices); 
in particular, 9 does not decrease. Since r = t\ + o(ti) as t\ — ► oo, we have proved 
the lower estimate in Lemma 1231 

2.18.3. A super solution. A slight modification of the preceding construction will 
give an upper barrier for q on the same domain f2. One begins with the a priori 
bound 

(56) u < C{1 + a 1 It) = (7(1 + p 2 ) and hence q < C{1 + p 2 )p 2 = C a (B)p 2 , 

which holds for A /^/t < p < B 2 . 

Choose a smooth function Q which satisfies Q(r) > r 2 /4, and define K n = 
max[ 0i i] \Q"\- Instead of requiring Q to be supported in [0, 1], we impose 

Q(l) = 2C (B) with C (B) as in 

For r > n, one introduces qo(p,r) = p 1 (r) 2 Q(p/p 1 (r)) , as in l|51|) . and finds that 
|-C[g ]| < K /t. For K > K and any (9 € (1, §), define 

q.+ (p, T ) = Q<ia{p, T ) - k/t. 

Choosing K > Ko + C(Aq,B), one finds as in (|53"|) that C[q+ — q] > for r > tq 
and Aj \fr < p < -B 2 . Thus 9+ is a supersolution, provided q+ > q on the parabolic 
boundary of f2. 

On the left edge of il, one has, as in (|55|l . g < (1 + Eq)A 2 /At. One also has 
g+ > (1 -ei)M 2 /4r, where ei = ei(A,if). 

Hence one will have q < q + on the left edge if 9 satisfies 

1 -±^<9< 3 -. 

1-ei 2 

On the right edge of one has p = Pi(t), and thus 

q+{p,T) = 9 Pi {t) 2 q(i) - k/r > 9(1 - ir/A 2 )g(i)p 1 (T) 2 g(iv 1 (r) 2 , 

provided A is chosen so large that 9(1— k/A 2 ) > 1. On the other hand we know that 
1 < C-o(B)pi(t) 2 holds on the right edge of O. Hence our condition Q(l) > Cq(B) 
implies that q + > q on the right edge. 

It now follows that q < q+ on f2, and as with the subsolution, this implies that 
q(p, t) < 9p 2 /A for all p G [A/^/t, B] and r > r x + o(n). 
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2.19. Convergence in the outer region. Let B and e > be given constants. 
Define A = A(e,B) and f = f(e,B) as in Lemma Then we have shown that 
when 

— B 2 1 

(57) a = By/r, i.e. when s 2 = —{T — t) log ■ 



4 v 7 b T-i' 
one has 

a/1 + (1 - e)B 2 /A <u< + (l + e)5 2 /4. 
For ?/> = R n (t)u, this therefore means that 



(58) y/T^ty/1 + (1 - e)B 2 /4 < - < y/T^ty/l + (1 + e)S 2 /4 

V 2 ( n ^ !) 

for all t € [i, T), where t = T -e^. 

We define the outer region (with parameter B) to be the portion of space-time 
given by 

r B = { (x, t) : f < r < T, s(x, tf > ^ (T - t) log ^ } . 

The boundary 8Tb of Tb is given by (|5T|) . 

For any given t* £ [i, T) , we let x* be the x coordinate of the point on dT b 
with t = t*. We will now show that s(x*,i) and i/j(x*,t) change very little for 
f* < t < T. This will directly imply the asymptotic description of VK 5 ? T) for s — > 
in the introduction. Fix some Bq > 1. We begin with a simple uniform estimate 
for ip s in the region B^y/r < a < rje ST . 

Lemma 24. For any Bq > 1, there exist constants 7, 6 > and C < 00 depending 
only on Bq and the initial data such that 

< M*,t) < ~F 

for all t > f(B ) and B ^/t < a(s) < ^e &T . 
Proof. Define r(t) = ip(0,t). By we have 



r= v ->(„ l) ( r-/i<jl-^ + (i) 



By we have 



> ^(«,t) > y/2(n-l)(T-t)yJl-B$/8 

for all s > s(Bq) — Bqt^JT — t and r > f(Bo) > tq. Now restricting our attention 
to the region of the neck where ip < 1/Bq and < ip s < I/-B0, we may apply 
estimate (34) of 0, obtaining C depending on Bq and F* such that 

logr 



#<Clog 



\ogip 



This estimate holds as long as a < je , i.e. as long as s < ji^/T — t) 1 < 
Let ci = log ^1 -B 2 /8. Then 



logr lo gv /2(7i-l)-§- jL+o(i) 2 Cl 1 

< . 21 — = H h O(-). 

logV' lo gv /2(n- 1) - § + ci r r 
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Hence r^ 2 < Cr log (gft) = C (2c, + o(l)) 



□ 



We are now ready to describe s(x*,t) and i)){x*,t). 

Lemma 25. For any x* the limits 

lim s(x#,t) and lim ib(x*,t) 
t-*T t->T 

exist. Furthermore, a constant C < oo which does not depend on B > Bo or x* 
exists for which one has 

(59) \s(x„t) -s(x„t,)\ < C{B- 3 + |log(T-t,)|- 3 / 2 } s(x„U) 
and 

C 

(60) \ip(x*,t) - tp{x*,U)\ < -^^{x*, t„) 



for allte [t*,T). 



Proof. To estimate the change in s(a;*, t) we use 
(61) 



ds 
dt 



J.{x*,t) | n f s(x " t] (tiY ds 



as in the proof of Lemma [21 Let t' e (i* , T) be given, and let x' be such that 
s(x',t') lies on 8Tb- 




(62) 



Then at time t' , Lemma implies that 

C 



< ^ s {s,t') < 



l0g ^7 



for all s < s' . 



Since x* > x 1 , we then get 



""ft)'*-, , 

V V / Js(x'A') V V 



^1 ds. 
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In the first term, we use to get 

^y*< - c r^^s 

ipj - y/\n\/{T-U) Jo V 
C 1 



/' 

Jo 



< 



< C 

" y/{T-t')\nl/(T-U)' 
In the second term, we have x' < x < x* so that 

i> > i>{x', t') > ^2{n - 1)(T - V) (1 + (1 - e)B 2 /4) > CBy/Y^V, 

with C independent of t or B > Bo- On the other hand, we also have ip s < 1, so 
that 

(z'.f) \V>/ ~ Js(x',t<) 

1 

< 



< 



c 



B^JT - t' 

For the first term in H61[). one has the same estimate, so that at any time t' 6 (t* , T] , 
one has 

< fta(**.0 < CiB- 1 + | log(T - i*)r 1/2 }(T - t')" 1/2 - 

This is integrable in t' G (i*,T), so that lim t ^T s(a;* , t) must exist. Integration 
with t± <t' <t then leads to 



\s(x»,t) - s(x*,U)\ < C{B- 1 + | log(T - U)\~ 1/2 }Vt^U. 

Recall apply a 2 6,a& 2 < a 3 + 6 3 to a = S" 1 and 6 = |log(T - t*)|" 1/2 , and 

conclude 

- s(x*,Q\ < C{B~i + | log(T - t*)r 1/2 }— ===== 

B^/lnl/(T - t*) 

< C{B- 3 + | log(T - <,)|- 3 / 2 } S (^, t„), 

which proves the estimate l|59|l for s. 

To establish the other estimate, we recall from Proposition 5.3 in |3] that there 
is a constant C which only depends on the initial data and the dimension n such 
that 

\2HH\ < C 

always holds. In particular, the constant C does not depend on our choice of B. 
In [3], we used this estimate to show that lim t ^T "^(x*, t) exists. It also implies 
that 

(63) \^{x^t) 2 -i){x*,U) 2 \ < C(t-t m ) 

and hence, using (|58p. that 
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for all t £ (t*,T). □ 

2.20. Asymptotics in the outer region. Using (|58|l and (|60|) . we compute that 
i){x*,t) = {l + 0{B- 2 ))i,{x,,U) 



= (1 + 0(B- 2 )) J2(n - 1) [1 + (1 + 9s)B 2 /4\ (T - U) 



for some \0\ < 1. 

From the definition H57J) . we see that s — s(x*, £*) satisfies 



By/T-t = , 

V-log(T-t) 

Hence, taking logarithms and abbreviating u> = T — t, 

log B + | In lu — log s — | In | log w| , 

so that 

logw = 21ns- log | log a; | - 21nB = O(logs) 
= 21ns + 0(log| log s|) 
= (2 + o(l)) logs 



We therefore find that 

BV^=(l + o(l)) 



V-21n; 

Thus as t* — > T, one has 

xl> = (l + 0(B- 2 + e ))(|Vn^T + o(l))- 



V-logs 

for all t £ [t*,T] and for which s = s(x* , t*) satisfies l|57|) . Since B and e can be 
chosen arbitrarily, this implies 

Lemma 26. For any 5 > £/iere exist _B > and i < T such that 



V-logs ^ V-logs 

in the region (T — f)| log(T — t)\ < s 2 /B 2 , t<t<T. 

3. Formal asymptotics for general neckpinches 

The usual starting point for obtaining matched asymptotic expansions near a 
stationary solution of a nonlinear pde is to linearize around that solution. This 
method encounters two obstacles when applied to the Ricci flow. 

The first challenge is the well-known fact that the Ricci flow is only weakly 
parabolic. This deficiency stems from its invariance under the full diffeomorphism 
group of the underlying manifold. We overcome it by a variant of the DeTurck trick 
[H1E3|- (See S E£H ) Our choice is equivalent to the Bianchi gauge, as adopted in 
the elliptic context by Biquard jHj and others. 
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The second challenge is more fundamental: it is the fact that the linearized 
pde possesses a null eigenvalue, suggesting the presence of a center manifold. This 
phenomenon also occurs in neck pinches of the mean curvature flow [S], and flat self- 
similar solutions to the reaction diffusion equation dtu — Au + u p |201 1211 l2~2l I13| . 
We overcome it by carrying out a higher-order variational analysis. A heuristic 
description of our method is as follows. We will formally compute the dynamics of 
a center manifold of a fixed point of a flow. That fixed point is the cylinder soliton 
introduced in § 13.11 and the flow is the dilated Ricci flow introduced in § 13.21 The 
center manifold will be tangential to the null cigcnspace of the linearization obtained 
in § 13.41 A general point on the center manifold will be of the form £ + $ (£) , where 
C belongs to the kernel of the linearization and <E> (£) is at least quadratic in its 
dependence on £. (See § 13.91 ) To compute the formal dynamics, we substitute this 
Ansatz into the dilated Ricci flow, project onto the kernel of the linearization, and 
compute the purely quadratic terms, yielding an ODE on the kernel. The errors 
introduced by this method should be smaller than the order of the solution, namely 
0(t~ 2 ), where r is the time scale introduced above in Section|21 

We pursue this method below. The computations it requires are extensive, but 
the (formal) conclusion it yields is eminently simple. It tells us that the asymp- 
totics obtained rigorously in Section [21 for rotationally-symmetric solutions should 
be stable for fully general solutions. 

3.1. The cylinder soliton. Recall that a Ricci soliton is a tuple (A4 m ,g,X), 
where M. m is a smooth manifold, g is its Riemannian metric, and A is a complete 
vector field on A4 m such that the identity 

-2Rc(.g) = £ x g + \g 

holds for some A e { — 1,0, 1}. It is well known that each Ricci soliton gives rise to 
a self-similar solution g (t) of the Ricci flow 

!*— 2Rc&) 

defined by 

g(t) = (T + \t) (fig), 
where ipt is the one-parameter family of diffcomorphisms generated by the vector 
fields (T + At) -1 A. 

Now consider the manifold RxS™ with local coordinates (y°; y 1 , . . . , y n ), where 
y° = x e R and (y 1 , . . . , y n ) Eye S n . Let g = g can denote the round metric of 
unit radius on S n . Define a product metric g on R x S n by 

(64) g = dx <£> dx + 2 (n - 1) g, 

noting that its Ricci curvature satisfies 

Rc(s) = (n-l)S. 

Let A be the vector field on R x S n defined by 

(05) _£|. 

noting that the Lie derivative of g with respect to A is 

Cxg = dx ® dx. 
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It is then easy to see that 

(66) -2Rc(g)=Cxg-9, 

hence that (R x S n ,g,X) is a shrinking gradient Ricci soliton. We call this the 
cylinder soliton. 

For later use, we observe that the Levi-Civita connection of g is given in local 
coordinates (y°; y 1 , . . . , y n ) by the Christoffel symbols 

(Pl7 \ r fc _ / i% ifl<i,j,k<n 

1 ' ^ \ otherwise. 

Here and throughout the remainder of this paper, a hat designates a geometric 
quantity associated to the round unit sphere (S n ,g). The components Rijke = 
gemRijk °f the Ricmannian curvature of g are 

(68) R = [ 2 ( n ~ 1 ) ^ii*< = 2 ( n _ l ) {9u9jk ~ 9ik9ji) i£ 1 < i,j,k,l < n 

3 \ otherwise, 

and those of its Ricci tensor are 

(69) Rij = 



Rij = (n — 1) gjk if 1 < i, j < n 
otherwise. 



3.2. The dilated Ricci flow. Motivated by the rigorous results obtained in Sec- 
tion [21 we want to study finite-time local singularities of the Ricci flow that are 
modeled on the cylinder soliton. Accordingly, let (M m , G (t j) be a solution of the 
Ricci flow 

^G = -2Rc(G) 
dt y ' 

that exists for < t < T < oo. Given any vector field X on A4 m , let ipt be the 
family of diffcomorphisms solving 

§i<p(;t) = T^-iX(<P (;*))• 

Then the blow-up of G at T modified by X is the metric g (t) defined by 
(70) G(t) = (T-t)(v*g(t)). 

Observing that 

| to) = T^Pl (Cx9) + & (j t9 ) , 

one computes that 

if* t (-2 Rc ( 5 )) = -2 Rc (G) = ^G = ^ (-g + C x g + (T - t) ^g 

Thus if one regards g as a function of the rescaled time variable 

1 



(71) r = log- r _ / . 
then (Ai m ,g (t) , X) becomes a solution of the dilated Ricci flow: 

(72) ^-g = -2Rc(g)~£ x g + g- 

The following is an immediate consequence of equations (|66|l and l|72|). 
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Lemma 27. The cylinder soliton (Ix S n ,g,X) is a stationary solution of the 
dilated Ricci flow. 

3.3. First and second order variation formulas. If ft is a symmetric (2,0)- 
tensor on a Riemannian manifold (A4 m ,g), its Lichnerowicz Laplacian A^ft is de- 
fined by 

(73) A e hij = Ahij + 2R ipqj hP q - R k t h k] - R*h ik , 
where Aft is the rough Laplacian. Its divergence Sh is defined by 

(74) (Sh), = -V'Ay. 

The following variation formulas are proved by direct calculation. 

Lemma 28. Let (Ai m ,g) be a Riemannian manifold; let ft be a symmetric (2,0)- 
tensor on M m ; and let e be small enough that 

g = g + eh 

is a Riemannian metric. 

(1) In local coordinates, the Christoffel symbols T of the Levi-Civita connection 
of g are 

y hi (Vi^ + Vjhi - V%) + O (e 3 ) . 

(2) In local coordinates, the Riemann tensor Rm of g is 
R l ijk = R ijk + |(ViV fc ft^ - V t V e h jk - VjVfeftf + VjV e h ik 

+ Y^(v 4 v m ft Jfc - VjV m h ik - v 4 v fe ft™ + v.Vfeftr) 
+ j 0W + v fc ft™ - v m ft 4fc ) (Vj/4 + v% m - v m /# 

- J (Vjft™ + V fc ft™ - V m ft, fe ) (V,ft £ m + V% m - V m ft^) 
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(3) In local coordinates, the Ricci tensor Rc of g is 

Aghij + ViVjH + Vi (Sh) j + Vj {6h) t 



—h pq {VpVqhij + ViVj-fcp, - ViVp/igj - VjVphqi) 



+ (h k R ]pkq + h k R ipkq + 2h k R ipk] ) + O (e 3 ) 

(4) //X and Y are vector fields on A4 m and X = X + sY , then 
C^g — £xg + e (Cyg + Cxh) + e 2 £yh; 

in particular, 



(L jt g) i . = V i X j + V j X i 



+ e [(ViYj + VjYi) + (X^kh^ + S7 l X k h k] + VjX k h ik )] 
+ e 2 (Y k V k h l3 + V t Y k h kj + V,Y k h lk ) . 

3.4. Arbitrary perturbations of a stationary solution. Let (A4 m ,g,X) be 
any stationary solution of the dilated Ricci flow l|72[l . and let h be an arbitrary 
symmetric (2, 0)-tensor such that the perturbation 

(75) g = g + h 

is a Riemannian metric on Ai m . Denote the trace of h with respect to g by 

(76) H = tx g h. 

Let Y be the vector field metrically dual to the 1-form \dH + Sh, namely 

(77) Y = (\dH + Sh 



2 

and form the perturbation 

(78) X = X + Y 

of the vector field X. The role of Y is to implement a DeTurck trick that 
makes the linearization of the dilated Ricci flow (|72[1 strictly parabolic. (In the 
elliptic context, this choice is frequently called the Bianchi gauge (5|.) Now applying 
Lemma |2H gives the following result. 

Lemma 29. The metric g (r) will be a solution of the dilated Ricci flow \7£fy for 
the vector field X if and only if h evolves by the nonlinear system 

(79) ^-h = Ah+Q(h)+C{h), 

OT 

where A is the linear elliptic operator defined by 

(80) Ah = A e h-C x h + h, 
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the quadratic term Q (h) is defined by 

(81a) [Q(h)l tj =h?V j V q h* + h$ViV«t% 

(81b) + h pq {ViV p h qj + VjVphqi ~ V p V q h i:j - ViVj/i M ) 

(81c) - - (h^VjVkH + tf\/iV k H + VitfVkH + V.,//fV/ //) 

(81d) + VitfVghl + VjtfVqhl + V p h q V q h p 

(81e) - Vphy p h jq - -V^f V^* 

(81f) + h pq [h^Rjpqk + hjR ipq k — 2h q Ripkj) > 



and C (h) is at least third-order in h and its covariant derivatives. 

The cubic term C (h) will not be important in the formal asymptotic analysis 
that follows. 



3.5. The linearization at the cylinder soliton. In the case that the stationary 
solution of the dilated Ricci flow is the cylinder soliton (Ix S n , g, X) introduced in 
the linear operator A defined in (|8U[) can be written in a more useful form. To 
see this, we again work in the coordinate system y , . . . , y n J , with y° = x E K 
and (y , . . . , y n ) = y E S n . Recalling formula (|69|l . it is easy to see that 

if i = j = 
if 1 < i < n 
if 1 < j < n 
if 1 < i, j < n. 

This observation allows one to write the Lie derivative of h with respect to X as 
Then using the simple identity 

RilRjk — RikRjl — — ^ — Rijkii 

one reaches the following conclusion. 

Lemma 30. On the cylinder soliton, the linear operator A defined in \8(]\) is given 
by 

(82) (Ah).. = Ahij - 1 + ^j R Pi ~ * Ai) 



!>', I'.i. ) — 




where Ah denotes the rough Laplacian. 
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3.6. An upper bound for the spectrum of the linearization. Since the man- 
ifold R x S n is noncompact, it is not necessarily the case that the formally elliptic 
operator A has discrete point spectrum. For our operator A the first order term 
— \xd x makes this happen, as we will now show. 

We consider the bundle 52 of (2, 0)-tensors on R x S n and the space C^°(6>2) 
of smooth compactly supported sections of this bundle. The operator A maps this 
space to itself. We define the inner product 

(h,k)= [ [ /iyfc tJ e^ 2/4 dO dx 

where d6 is the volume form on S n . We write \\h\\ = \/{h, h) for the corresponding 
norm. The completion of Cq 30 ^) with respect to this norm is L 2 (S2', e~ x / 4 d6 dx), 
which we abbreviate to L 2 . (See for instance.) 

Lemma 31. The closure of the densely defined operator A is a self-adjoint operator 
with compact resolvent. A is bounded from above. Its spectrum consists of discrete 
finite multiplicity point spectrum. 

Proof. One can write A as 

(83) (Ah) tj = e * 2 / 4 |_ { e -* 2 A} + ^Ah l3 + V ijpq h™ 

in which A is the rough Laplacian in the direction tangential to the spheres {x} x S n 
and 

Vijpq = ■^ZT[(RijRpq ~ RiqRjp)- 

By integration by parts, one finds that for any h,k € Cg ^) one has 

(84) (h,-Ak) = 

{O x hijd x kij -\- 2(n—l) ' V gk{j Vijpqh{jkpq} 6 ^ dO dx, 

is™ 

from which it is evident that A is symmetric. Boundcdness of Vij pq implies that 
there is a constant Cq such that 

(85) {h,~Ah)>J J \Vh\ 2 e- x2 ' 4 d6 dx~C \\h\\ 2 

holds for all heC^°(S 2 ). 

Let H 1 = H 1 (S2) be the completion of C^°(6>2) for the norm corresponding to 
the inner product 

(h,k) = / / {Vih pq Vik pq + h pq k pq )e~ x /4 dO dx. 
Jr Js n 

(Again, see [H].) Then the standard Hilbert space arguments show that for all 
A > Cq, the equation Xh — Ah — k has a distributional solution h G H 1 for any 
k € L 2 . It follows that the closure of A is a sclf-adjoint operator in L 2 whose 
domain 

D(^) = {(A-^)- 1 fc|fcei 2 } 

is contained in H 1 . (Here it makes no difference which A > Co one chooses.) In 
fact, the inequality l|85(l implies that the H 1 norm is equivalent with the norm 
2 - (h, (A -A)h), so that H 1 is the form domain of A, i.e. H 1 = D(VA - A) 



SIGURD ANGENENT & DAN KNOPF 



The inequality lj%5|) tells us that the spectrum of A is contained in the interval 
(-00, Co]. 

To see that A has compact resolvent, we recall Lemma which implies 

\\xh\\ 2 <cf f {\d x h\ 2 + \h\ 2 }e~ x2/4 de dx <C\\h\\ 2 Hl . 
Jr Js" 

Hence 

\\h\\ 2 + \\xh\\ 2 + \\Vh\\ 2 < C(h, (A - A)h). 
This inequality together with the Rcllich-Kondrachov theorem imply that the imbed- 
ding H 1 C L 2 is compact. Since H 1 is the form domain of A, we conclude that A 
indeed has compact resolvent, so that its spectrum is pure point spectrum. 

□ 

Remark. More generally, the linearization of the modified Ricci flow operator 
at any gradient Ricci soliton (A4 n , g) will be self-adjoint in L 2 (A4 n , g; d/i) if dpi 
is the measure naturally associated to g and one takes ip to be the soliton potential 
function. 

3.7. A decomposition into invariant subspaces. In order to compute the spec- 
trum of A, we split the space L 2 (5 2 ) into a number of A- invariant subspaces. 

It is a standard fact that any (2, 0)-tensor h on a Riemannian manifold (A4 m , g) 
admits the decomposition 

h = h v +h°, 

where 

/i v = {tr g h) g 

is a multiple of the metric, and 

1 



h° = h-—h v 

777 



is trace free. 



One can exploit the product structure of the cylinder soliton (K x S n ,g,X) to 
obtain a more refined decomposition. Again using coordinates (y°; y 1 , . . . , y n ) and 
adopting the convention that Roman indices lie in the range . . . n while Greek 
indices lie in 1 ... 77, we decompose a given symmetric (2, 0)-tensor h on K x S n as 

(86a) h = M dx®dx + N a (dx ® dy a + dy a ® dx) + (O^ + 0° a/3 ) dy a ® dy 13 

(86b) = M dx®dx + N a (dx <g> dy a + dy a <g> dx) + (Pg a/ 3 + Q a p) dy a ® dy . 

Here one should regard M as a scalar function on Kx S n , N (x, •) as a one-parameter 
family of 1-forms defined on the spheres {x} x S", P as the scalar function 

77 

and Q (x, •) as the one-parameter family of trace-free (2, 0)-tensors defined on the 
spheres {x} x S n by 

Qa/3 = h aj 3 — Pcja/3- 

The utility of this decomposition is that it exhibits a set of invariant subspaces 
for A. To demonstrate this, it will be helpful to fix additional notation. Let A 
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denote the rough Laplacian of the canonical sphere metric g. Let —AdR denote its 
Hodge-de Rham Laplacian, recalling that AdR acts on a 1-form N by 

A dR N a = - [{dS + Sd) N] a = AN a - RPNp = AN a - (n - 1) N a . 

Let Ai denote the Lichnerowicz Laplacian of <?, which acts on a trace-free tensor Q 

by 

AeQ a /3 — AQ a p — 2nQ a p. 
Finally, let B denote the differential operator defined by 



(87) 



B = 



if 











d 



dx 2 2 dx dx" dx' 

Using these conventions, one computes using equation (|82|) that the decomposition 
of Ah corresponding to Ip3|l is 

1 



C^JOoo = 



B 



2(n- 1) 
1 

2 (n - 1) 
1 



2 (n - 1) 
1 



2 (n - 1) 



A 



M 

N a = 

A + l 
1 



B 



n- 1 



Qa/3- 



This calculation proves the following result. 

Lemma 32. On the cylinder soliton, the operator A given by \8ty) may be repre- 
sented schematically with respect to the decomposition 1^86)) as 



A 



(M\ 

N 
P 



I 



B 



4- 1 

T 2(n-l) 



A 



M 



\ 



B 



2 („-l) A rffl + 2 



N 



B 



2(n-l) 



A + l 



p 



B 



i 

n-l 



QJ 



2(ra-l) ' 

The only contribution to the spectrum of A comes from its four component opera- 
tors: B+ 2 ( n 1 _i) A and B+ ^n-i) A + l acting on scalar functions, B+ 2 (n-i) ^-dR+\ 
acting on 1-forms, and B + 2 (n-i) ^ ~ TT^T ac ^ n <7 on trace-free (2, 0)-tensors. 

3.8. The spectrum of the linearization. The decomposition obtained in Lemma 
[fallows us to analyze the spectrum of A using separation of variables. To illustrate 
the idea, suppose that 

M {x, y) = h (x) F (y) , 
where h is an eigenfunction of B with eigenvalue A, and F is an eigenfunction of 
2(n-i) w ith eigenvalue /j,. Then 

1 



B 



2 (n-l) 



-A 



M = (A + (J,) M. 
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This line of argument leads easily to the following observation. 

Lemma 33. All eigenvalues of A : D(„4) — > L 2 have the form A + fi, where A is 
an eigenvalue of B acting on f) = L 2 (M.;e~ x / 4 ), and /i is an eigenvalue of either 
2 (n-i) A or 2(n-i) ^ + 1 actin 9 on C°° °f 2 (n-i) Arffl + \ acting on fix (S n ), 

or of 2( w 1 _i) A ~ ^rr acting on the space of smooth trace-free (2,0)-tensors on S n . 

To make this result useful, we need some simple facts about the component 
operators that appear in l|88(l . 

The operator B = — § ^ acting on scalar functions. The Hermite 
polynomials {/ife}^ constitute a complete orthogonal family for the operator B 
on the weighted Hilbert space 

f) = L 2 (R,e~ x2/4 dx^ . 

Moreover, one has 

(B + l)h k = h k , 

so that the spectrum of B is {0, — h, — 1, — §, . . .}. (Notice that f) is a larger space 
than we defined in Section[21 since we not restricting to even functions here.) 

The operator „/ 1 lS A acting on scalar functions. It is well known that 

2(n— 1) 

the spectrum of A acting on C°° (S n ) is {Afc} fc>0 , where 

^fe = — k (n + k — 1) . 

In particular, the only nonnegative eigenvalue of 2 (n-i) * s A ° ~ 0- ^ s eigenspace 
consists of the constant functions. 

The operator 2 ( ra 1 _ 1 ) A + 1 acting on scalar functions. If k > 2, then 

A fc + 2 (n - 1) = (2 - fc) (n - 1) - k 2 < -k 2 < 0. 

Hence the only possible nonnegative eigenvalues of 2 (n-i) ^ + arc 

A + 2(n- 1) 
2(n-l) 

and 

Ai +2(n- 1) _ l n-2 f 1 
2(n- 1) ~~ 2n - 1 G [ ' 2 
The eigenspace corresponding to 1 consists of constants. The eigenspace corre- 
sponding to i ^5? consists of the spherical harmonics: the restrictions to S n of the 
linear functions on 

The operator 2 ^_^ A d R + \ acting on 1-forms. It is shown in [7] that 
the spectrum of K^r acting on the space Qi (S n ) of smooth 1-forms is 

{^} fc >iUi Afc + 2 -"}fe>i- 

If k > 1, then Afc + n— 1<— 1. Because Afc + 2 — n < A& in all dimensions n > 2, 
it follows that 2 (n-i) ^dR. + 5 is negative definite. 
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The operator „, 1 n- A acting on trace-free (2,0) -tensors. The 

2(n— 1) 71—1 v ' ' 

operator 2 (n-i) ^ ~~ T^i ^ s cl eai "ly negative definite. 

Combining these observations with Lemma 123 leads to the following conclusion. 

Lemma 34. On the cylinder soliton, the only possible non-negative eigenvalues of 
the linearization A are 0, and 1. 

7 2 n— 1 7 2 7 

• The eigenspace corresponding to 1 consists of constant multiples of g. 

• The eigenspace corresponding to h consists of multiples of g that are linear 
in x G R. 

• The eigenspace corresponding to i consists of spherical harmonics. 

• The eigenspace corresponding to consists of constant functions and of 
multiples of g that are quadratic in x £ R. (When n = 2, the spherical 
harmonics also belong to this eigenspace.) 



Remark. The results above can be obtained in another way. Since all geometric 
data of the metric l|64|l are independent of x £ R, one has the simple commutator 

' d] _ ld_ 
' dx\ 2 dx 

So if Ah — Xh, a straightforward induction argument shows that 



dx k J \ 2 J \dx k 

for all integers k > 0. Since Lemma 1311 proves that the spectrum of A is bounded 
from above, it follows readily that for every eigentensor h of A, there exist some 
integer k = k (h) > and a family of (2, 0)-tensors hi, . . . ,h K independent of x E R 
such that 

K 

h = 22 xk hk- 

k=0 



3.9. Our Ansatz. The analysis in § 13.81 shows that any eigenvalues of A corre- 
sponding to trace-free (i.e. non-rotationally symmetric) perturbations of the cylin- 
der soliton are all strictly negative. As we explained in the introduction to this 
section, the only perturbations relevant to our formal center-manifold computation 
are those corresponding to null eigenvalues of the linearization. Therefore, we adopt 
the Ansatz that 

h = M dx <g> dx + N a (dx dy a + dy a ® dx) + (Pg a /3 + Q a p) dy a (g) dy 13 
has the special form 

(89) h (x, y,r) =u (x, r) dx ® dx + {v (x. t) + w (y, r)} g, 

where once again, r denotes rescaled time, x denotes a point in R, and y denotes a 
point in 5*™. Specifically — continuing to let {hk}^L denote the Hermite polyno- 
mials — we assume that M = u, where 

(90) u (x, t) = u (t) h (x) = u (r) . 

We further assume that N = 0, that Q = 0, and that P = v + w, where 

(91) v (x, t) = v 2 (r) h 2 (x) = v 2 (r) (x 2 - 2) 
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and 

(92) w(v,t) = { Si=i"iW ni(y) ifn = 2 
v ; 7 \ otherwise. 

Here (f2i, , SI3) is a basis for the space of spherical harmonics on S 2 . (Recall 
that w corresponds to a positive eigenvalue of the linearization whenever the total 
dimension n + 1 is at least 4.) 

To implement our formal center-manifold analysis, we shall study the flow 

(93) ^-h = Ah+Q(h) 

or 

that models equation (|79|l up to second order. To do so requires us to analyze 
the quadratic term Q (h) at the cylinder soliton in a manner analogous to what 
was done above for the linearization A. To simplify the notation, we again use 
coordinates (y°; y 1 , . . . , y n ) and assume that Roman indices lie in the range . . . n 
while Greek indices lie in 1 . . . n. If / is a smooth function, we further adopt the 
convention that f x = df /dx and f a = df /dy a . 

We shall begin with computations that are more general than what we need at 
the moment. To wit, we assume only that h has the form given by l|89|) . without 
imposing the specific assumptions 19U|) - (|92|I . One then verifies readily that all first 
covariant derivatives of the tensor h vanish except 

Vo^oo = u x 
Vo/la/3 = v x g a @ 
V~(h a p — w~f g a p. 

All second covariant derivatives of h vanish except 

VoVo/loO = u xx 

VoV /i a /3 = v xx g a/ 3 
Moreover, the only nonvanishing derivatives of the trace 

■ Ti 

H = g ll h l3 =u + - _ (v + w) 

are 

n. 



2(n-l) 
n 



7 ~ 2(n-l) 7 
and 

V V i? = Us 



2(n-l) 



n 



V 7 V 4j ff - — -V 7 V 5 u>. 

2 (n — 1) 
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Using these formulas, a straightforward calculation reveals that the first compo- 
nent of the quadratic term 

[Q (h)] 00 =2/i 00 VoV /ioo - h at3 V Voh a f3 - h 00 VoV H 
- Voh Q oVoH + | (Vo/ioo) 2 - ivo^Vo^ 

may be written as 

(94a) [Q (h)] 00 =uu xx -^v xx [2(n -l)u + v + w] 

4 (n — 1) 

(94b) + X -ul - —?— v x [4(n - 1)^ + v x ] . 

Then using the fact that 

g^ 5 %V s w = 1 Aw = - - -r ti>, 
2 (n — 1) 2 (n — 1) 

one computes that the remaining components 

[2 (h)] aP =hZVpV s h s 7 + ^V Q V*/i 7 - h 00 V \7 h af3 
+ h^ s (V Q V 7 

+ V a hZV 5 ht + VphlVsh* + Vjh 5 a V s h} 



h lS (hlR 0lSv + h}R 



are given by 




(95a) [CWU=- 


1 


2(n- 1) 


(95b) 


1 


4(n- l) 2 


(95c) 


n-2 
8(n- l) 2 


3.10. The inner layer |x| 





Vw 



v + w) \7 a \7 pw + 3WaWp\ 



9 {x, V,r)=g (x, y) + h (x, y, r) 

near a singularity modeled on the cylinder soliton g. We assume that (18911 and 190(1 - 
(I92[l hold. These assumptions force h to belong to the kernel of the linearization of 
equation lf7U|l . because 

Ah = [Bu] dx®dx+ [{B + 1) v] g = 0. 

Recall that w vanishes by assumption if n > 2. In this context, Q94J) becomes 

(96) [Q (h)] Q0 = - \ v 2 w + - 2 [(n - 1) u t> 2 + vf] L - J - 2 v| I fo. 

And (|95|l reduces to 

[Q(h)U = Pg a p, 
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1 












4 





2 4 1M 2 

2 I u i J „,. „.2 



where 

- 2u w2 —r v 2 } h o + { 7;v 2 w — ^ /i 2 . 

§ n — 1 ) [2 n — 1 J 

In order efficiently to derive ODE for the functions uq, v 2 , and wi, u/2,ct>3 from 
equations l|96|) and l|97|l . it is helpful to adopt some additional notation: if / is a 
function of y € S n , we define 



ll/lls" = / f(y) dm 

and if X and Y are symmetric (2, 0)-tensors 01 R x S", we define 

/>oo 

((X,Y)) = / (X,Y) x e-* 2 / 4 dx, 



where (X,Y) X = fr x y xS „ (X,Y)~ dp,. By using the fact that {hk}^ =0 is an orthog- 
onal family for f), it is then easy to see that 



d 



«o(rHIM5l|i|ls« = Jt ((h,h dx 2 )) 

= ((Q(h),h dx 2 )) 



(n-1) 

and 



2 [(n- l)u ^2 + ^] ' ||J|ls- 



v' 2 (r)-n\\h 2 \\l\\lf sn = ^ t ((h,h 2 g)) 

= ((Q(h),h 2 g)) 



^ 2 11. 1 1 2 in ,|2 

"«2 • n \\ h 2\\ 



Similarly, for any i € {1,2,3}, we find when n = 2 that 

<4 ( T )- 2 \\ h o\\l \\£U\\ 2 s» = j t {h,siig)) = ({Q(h),n i9 }) =0. 

This calculations let us derive a formal profile of h in the inner layer |x| = o (y/r). 
Lemma 35. For perturbations of the form 

h = {u (r) h (x)} dx ® dx + |w 2 (r) /i 2 (z) + (r) f2j 5, 

£/ie /Zow * s equivalent to the system of ODE 

(98) -r M o = - ~ . 2 [O ~ 1) ^2 + v 2 ] 
ar (n — 1) 

(99) ±V2 = - —vt 
dr ri—l 

(100) ^ Wl = (* = 1.2,3) 
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whose solutions, up to an error term of o(t 1 ) as t — > oo, are 
(101) uo(t) = 



2t 



ra 1 



otherwise 



2(n-2) t 

(102) r, 2 (r) = 2^1 

(103) w i (r)=w i (0) (i= 1,2,3). 
Formally, this is the profile of h in the region \x\ — o(y/r). 



Proof. If «2 is not identically zero, equation (|99|l has the explicit solution 

n — 1 



u 2 (r) 



c + 2t 

Then equation becomes 

d ( n \ n 

— u Q It) + — u = 77. 

dr \c+2t) (c + 2t) 2 

An integrating factor for this linear ode is (c+ 2r)™^ 2 , whence one obtains 



"o (r) 



. log(c+2r) g if 2 

c+2r ' c+2r 



-^—r^ — I- -, — C otherwise. 

1-2 c+2t (c+2t)"' 2 



a 



Remark. To compare the formal results above with those obtained rigorously 
in Section |5J observe that our choices of dilating factors imply that the quantity 
a 2 (t) appearing in Section [21 and the quantity v 2 (r) appearing above are related 
by 

I \ ~ ^ (t) 1 

4 [n — 1) ot 
because in the inner layer, we have x w a and 

{2{n-l) + v 2 (r) ^ 2 (o;)+o(t- 1 )} .g 

= 2 (n - 1) [1 + a 2 (t) ftj, (a) + o^ 1 )] 2 g 

= {2 (n - 1) + 4 (n - 1) a 2 (r) fc a (<r) + o^ 1 )} 5. 



Remark. By modifying g by an initial conformal diffcomorphism, one may 
make all uji vanish. Therefore, throughout the remainder of this paper, we will 
assume that u>i = for i — 1, 2, 3. 
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3.11. The intermediate layer |x| = O (y/r). In order to study the region where 
|ir| = O (y/r), we replace the coordinates (x, t) with (£, r), where 



(104) 



If / is a smooth function of (x, r), we define F (£, r) = / (x, r) . Then one has the 
formulas 

e _ i 



fr=F T 



2t 



-Ft 



and 



fx 



--F f 



Let 



U (£, r) = u (x, t) and V (£, r) = u (x, r) , 
where u and w are the quantities that appear in the simple Ansatz i|89|) . Then 
substituting formulas l|94l) and (|95|l into the PDE l|93|) . one computes that 

(i + to % + 1^ + \u 2 

U T = -^U 6 + -\ -jr^V^{n-l)U + V) 

Z T 



^^[4(n-l)t/ e + y,] j 



and 



i 



V T = V - j-Vj: + - < V a (1-U) + ^ 



1 yj 
2 (n - 1) 4 



Recalling Lemma 13*51 we will construct an approximate solution by considering the 
equations 

L 



and 



V t = V-^V 6 + 0(t~ 2 ). 



These equations suggest that the solution is modeled by the time-independent pro- 
files U (£) = C\ and V (£) — C^ 2 , respectively. Matching this layer with the results 
of Lemma 1351 determines C\ and C2, yielding the following result. 

Lemma 36. In the intermediate layer \xj yfr\ = |£| = O {1), the formal solution is 
modeled by 

n — 1 x 2 



u(x,t)mJ7(-^) =0 



v(x,T)f*V[-= = 



To wit. 



g « dx <g> dx + 2 (n - 1) 1 + (£/2) 2 



3.12. The outer layer |a;| = o (l/\/T — £)• To describe the outer layer, we first 
introduce the "blown-down" coordinate 



(105) s = e~ T/2 x = y/T-t ■ x, 

which should be compared to the metric distance s defined in Section [21 

If one were to follow the method above by defining U (s, r) = u (x, r) and 
V (s, t) — v (x, t), an easy computation would show that 

U T =G{T-t) and % = V + O (T - t) . 
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Unfortunately, the steady-state solutions to these equations fail to convey enough 
useful information. 

Instead, we proceed as follows. Fix a large number 5 > 0. For each s, let t s be 
the time such that |£| = H at t s , where £ is defined by l)l()4[l . Let r s = log jhr- 
Then 

3= i^i = _w_ ,s| 



fr. ^(T~t s )r s 
This shows that 

T a = log (S 2 ) + log — + logr s , 

hence that 

(106) t s = [2 + o(l)] log ^ as |s|-»0. 

Now we "blow down" the solution of the dilated ( "blown- up" ) Ricci flow by defining 

u (s, r) = e~ T u (x, t) and v (s, r) = e~ T v (x, t) . 
Then at time t s , we compute that 

, ■'■ \ n-ls 2 
v[s,Ts) = e "V 



Ta J 2 r. 



When |s| is small and S is large, the undilated solution v cannot change much in 
the short time t s < t < T . Hence (jlOfill implies that as \s\ — > 0, one has 

V(S,T S ) = 3 rl± + 0(l)] 

2 [2 + o(l)] log 
n — 1 s 2 . , 
2 log^ 

Since u (s, t) = er T U (x/y/r) = 0, we have obtained the following result. 
Lemma 37. In the outer layer \x\ — o f -^== ^, the formal solution is modeled by 

u{x,t) w ^-u(«,t) anrf u (x, t) w v (s, r) , 

w/iere 

— 1 s 2 

u(s,r) = and v(s,t) = =— [1 + o(l)l . 

2 log 

To wit, 

. . n — 1 x 2 
u (a;, r) 



2 r + log ■ 
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